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We show that there are implications on thermodynamics 
that come from the existence of the initial cosmic singularity. 
At present time this is more a conceptual change than an ob- 
servable one. However at very early cosmic times there is a 
big difference between the actual behavior of thermodynamic 
quantities and the behavior assumed in the standard cosmo- 
logical model. We present the discussion of two systems: an 
ideal monatomic gas at present, and a photon gas at the early 
Universe. We show the striking result that the entropy den- 
sity goes to zero as the cosmic time goes to zero. This in turn, 
provides an explanation for the so called horizon problem. 

PACS number(s): 98.80. C, 05.20, 65.50 



When one studies a thermodynamic system one nor- 
mally have the thermodynamic laws, resting on our 
knowledge of the phenomenological equations of state; 
and one also seeks for an explanation of all these laws 
from statistical mechanics. If for example we study a sim- 
ple system like a classical ideal monatomic gas com- 
posed of N molecules in a volume V at temperature T; 
we have the phenomenological thermodynamic descrip- 
tion in terms of the equations of state for the pressure 
P and the internal energy U ; namely Pph = {NksT) /V 
and Uph — ^NksT [|l|; where ks is Boltzmann's con- 
stant and we are using the subscript ph to emphasize that 
these expressions have a phenomenological origin. One 
also has the explanation of these equations from a mi- 
croscopical description from statistical mechanics. In the 
process of performing a calculation in the framework of 
statistical mechanics, one needs to settle the appropriate 
distribution of probabilities for the microscopic states. 
This distribution is normally calculated from the maxi- 
mum entropy principle Q ; which basically states that 
the probabilities pi for the system to be in the micro- 
scopic state i is to be calculated by maximizing the en- 
tropy S{U, V, N) subject to the information that we have 
on the system. For example in the above case of the clas- 
sical ideal gas, one would say that while the volume V 
and the number of molecules N are fixed; we only know 
for the microscopic energy that its average value must 



coincide with U; that is, J^Pi^i = '^'i since the fact that 
our system is at temperature T, precludes us from a more 
precise information on the energy of the system. There- 
fore to calculate the distribution of probabilities {pi} we 
must maximize the entropy 



5* = -fc^y^Pt In Pi 
subject to the condition 

J2P^E, = U. 

The distribution so obtained is given by: 



(1) 



-PEi 



-/3iS, 



(2) 



where /3 = l/^ksT); which is recognized as the Gibbs or 
canonical distribution. 

We have said before that the entropy must be max- 
imized subject to the information that we have on the 
system; and we do want to take into account that we live 
in a Universe with an initial cosmic singularity. The way 
in which this observation affects our statistical mechanic 
calculation is by imposing some constrains on the avail- 
able phase space of our system. More precisely, the ex- 
istence of an initial cosmic singularity implies that there 
are particle horizons |^ ; which means that our past and 
the amount of matter that we can observe are finite. 
The implication of this is that we can not assume that 
our thermodynamic system is in contact with an infinite 
reservoir at temperature T. In particular this means that 
not all of the theoretically predicted energy levels Ej are 
available; since, if we label with -Emax the total avail- 
able energy in our causal past, then clearly we can only 
consider Ej < E^i^^. 

Following the standard calculation for a classical ideal 
gas we find that the partition function Z — J^j e~^^' , 
appearing in the denominator of (|^) above, is given by 
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where is the energy of the particle i in the state j. 
In the usual presentation of a classical ideal gas, it is 
assumed non-interaction and statistical independence for 
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each atom of the gas; from which it is deduced that the 
sums over each of the indices appear as factors; so that 
the partition function would be given hy Z = /Nl, 
where z is given, in the continuous representation, by 
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where h is Planck constant. The correcting factor 
accounts for the indistinguishability of the particles |^] 

Due to the restriction on the phase space coming from 
our cosmological knowledge we are not free to integrate in 
up to arbitrary large values of p^; since ^ can not be 
larger that the total available energy i?max m our causal 
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past; in other words we must have < -Emax- There- 



fore, one has z — ^ 
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"^x^/'^dx. In 

the limit for -Emax going to infinity one obtains the usual 
result Zoo = {V/h'^){27rmkBT)^^'^; however due to the 
cosmological restriction we can show that 

\nZ{(3,V,N) = N\n(^^^2^mf^ ~^N\nf3 

+Nh). / e-''x^/'^dx-\nN\ . 
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Then the entropy is given by 

S ^ kBl3^—^E, + kB \tlZ ^kePU + kB \nZ. 
And the corresponding energy is given by 



U = NksT 
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with Xc = /3£'max- To givc an idea of the order of 
magnitude of the correcting term, let us consider T = 
300°K, and a cosmic mass density p to be two per- 
cent of the critical mass density; that is, p = 0.02pc = 
(0.02)3i7^/87rG, with the Hubble parameter with 
value Hq = 75{km/s)/Mpc; and where G is the gravita- 
tional constant. Assuming a sphere of radius R = c/Hq, 
where c is the velocity of light, one can estimate -Emax 
by E-[aa.-K — 4,TT pc^ / 3; which implies a value of Xc on 
the order of Xc — 10^*. Such huge value for Xc gives a 
negligible correction in front of 3/2 in equation Even 
if one considers temperatures as high as those found in 
the center of the sun, one will have x^, ~ 10^^; and the 
whole correcting term to the equation of state , would 
be of order 10~^^^ in front of the 3/2. 
The pressure can be calculated from 

IdlnZ _ TV 
^ P dV ~ 

of state for the pressure 



therefore in this case the equation of state for 1 
does not show any correction; that is P = Pj 



It is observed that in the previous expressions the cor- 
recting term is more important if /^-Bmax becomes smaller. 
This situation arises when one approaches the initial cos- 
mic singularity; since in this regime /3 goes to zero and 
also i^max shows the same behavior; due to the fact that 
the particle horizon reduces its size as one considers ear- 
lier cosmic times. 

At very early cosmic times, during the radiation dom- 
inated era, the different kind of particles contribute to 
the energy momentum tensor as different components of 
a relativistic gas. In order to be concrete we will study 
next the system corresponding to a photon gas in a 
small volume V at temperature T, at very early cosmic 
times. 

The corresponding phenomenological equations of 
state for this system are: Uph ~ a V T'^ and Pph = 
\Uph/V\ where a ~ '^crvi/c, and cto is the Stefan- 
Boltzmann constant. From the knowledge of the equa- 
tions of state, one can calculate the phenomenological 
entropy; which is given ||] by Spu = 

In a Friedman Universe the line element is given in 
terms of the metric of a Robertson- Walker spacetime, 
namely H| 



ds^ = dr - 
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dl 
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where A(t) is the expansion parameter with units of 
length and dl'j. is the line element of a homogeneous and 
isotropic space with constant curvature fc = 1, or — 1. In 
a radiation dominated Universe, pA^ is a constant; and 
for very early cosmic times, A{t) behaves asymptotically 
as -v/t; where the initial singularity is attainned in the 
limit t ^ Q. Then, one can deduce that the phenomeno- 
logical entropy density, sph = Sphiyo)/Vo, diverges as 
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Let us now proceed with the correct causal statistical 
mechanic calculation of the entropy, taking into account 
our knowledge of the initial singularity. As we have seen 
in the previous example, we expect to find a different 
expression for the correct entropy; which implies that 
the corresponding equations of state will be different to 
those that appear in the standard cosmological model 
dominated by radiation. This means that to be abso- 
lutly consistent one should solve the Einstein equations 
taking this effect into account. We plan to do this in a 
future paper; instead here we would like to concentrate 
in the general behavior of the entropy in a Universe with 
an initial cosmic singularity, by given the calculation in a 
fixed geometry. We choose this geometry to be the Fried- 
man line element corresponding to a Universe dominated 
by radiation. We also have the freedom to relate one of 
our thermodynamic variables to the geometry; we do this 
by demanding the energy density to behave as l/A**; as 
it does in the standard cosmological model of a Universe 
dominated by radiation. 
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In order to be explicit, let us consider an open Uni- 
verse; that is, the k = —1 case, whose expansion param- 
eter is deterimined from the Friedmann equations H 




(7) 



and where the density p is linked to the geometry in 
such a way that the cuantity Qr, defined by Qr = 
{8ttG/3)pA'^, is a constant. 

Let us consider a reference volume Vq, which is small 
enough so that any other dynamical time is much larger 
than the typical thermodynamic relaxation time |^] asso- 
ciated to our subsystem determined by Vq. An observer 
that is comoving with the matter flow ascribes the en- 
tropy S{t) — S{U{t), Vq) to this small subsystem; where 
the world line of the observer can be parametrized by the 
cosmic time t. Let us recall that the chemical potential 
for a photon gas is zero; therefore the entropy does not 
depend on the number of particles. 

Following the standard procedures appearing in text 
books ||] of statistical mechanics one can prove that the 
entropy and internal energy for electromagnetic radiation 
contained in a volume Vq are given, in the continuum 
representation, by 



Vnk 



OKB 



4 

+ 3 









dw 





w 

T 



and 



U 



Vo 



{e^ - 1) 



dw : 



where a;„iax(i) = Pit) Emax{t) and the dimensionless vari- 
able of integration w is defined hy w = phw. 

If there where no particle horizons the first term would 
not contribute to the value of the entropy; since in this 
Then, in this situation one would repro- 
duce the expression for Sph mentioned above. 

It is important to remark that the Stefan-Boltzmann 
law; which gives the energy radiated per second per 
unit area / for black body radiation, namely, / = 
{cUph)/{4:V) = (To T'^, is now corrected by the expres- 
sion (|9|). One way to understand this result is that the 
Stefan-Boltzmann constant ctq, mentioned above, is ac- 
tually a parameter that changes with cosmic time; which 



is given by a{t) 



I 



(e--l) 



dw; and where 



one usually approximates for the present time to by the 

j,4 °o r 3-1 

value cr(to) ^ ao = , og^^ / j-^-pr dw 



When constructing a cosmological model, one normaly 
represent matter by some kind of fluid. In any Uni- 
verse whose energy momentum tensor is that correspond- 
ing to a fluid, given an event p, one can calculate the 
total available energy in the causal past of p by inte- 
grating the contribution of the fluid on the boundary 
of the causal past of p. In our case we can calculate 
Emaxit) explicitly by considering the contributions on 
the past light cone at earlier times t' coming from the 
energy density p{t')c^, with the corresponding red shift 
factor A{t')/A{t) and multiplied by the volume element 
47r r{t,t')^A{t') dl = Att r{t,t'fc dt'; where r{t,t') is the 
radius of a sphere, with center at the point of observa- 
tion at cosmic time t, and with surface on the past light 
cone at time t' . Considering the realistic case fc = — 1, 
corresponding to a Universe with energy density below 
the critical value, the line element of the three dimen- 
sional homogeneous and isotropic space can be expressed 
by dP = dx^ + sinh(x)^(d6'2 -|- sm{9)'^ d(p'^) . Using these 
coordinates, the radius r{t,t') of the sphere is given by 
r{t,t') — A{t') smh{x{t,t')); where the center of the co- 
ordinate system is taken at the point of observation, so 
that xit^t) = 0- the spacetime diagram in Fig. |l] one 
can see a picture of the coordinate system. 



observation point 



direction (9, <|)) 




past light cone 



initial singularity 
FIG. 1. The observation point is the tip of the past light 
cone with the generic time value t. A direction of observation 
with angular coordinates {8, 4>) reaches the sphere in the past 
at time t' < t, with radial coordinate X- The past light cone 
extends up to the initial cosmic singularity with the time value 
t = 0. 

One can calculate x{t, t') from the expression x{t, t') = 
c Jj, dt" /A{t"); from which it is obtained Ema.x{t), namely 

i?max(t) = f p{t'y^4nr{t,t')^c dt' 
Jo ^(t) 

3QrC^ 1 



(t+ 1) + 27 



iG A{t) 
- ln(r + 1 + Vt^ + 2r) 



where we are using the convenient dimensionless time 
variable t = {t(? j \fQ^). In order to study the behavior of 
Emax{t) near the initial singularity it is useful to express 
it in terms of A{t)] which is also a good time variable, 
since it is a monotonically increasing function of time. In 
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this way one obtains 




It is clear that as t decreases, so does the region covered 
by the particle horizon; and therefore E„^^^{t) also de- 
creases. More specifically one can show that the asymp- 
totic expression of ii^max for small values of A starts 
with: E^^^{A) = (3\/Q;c3/4G)[(2/3)(c^/VQ7)' + 
0{{cA/ ^fO^)"^)]. In particular, one can see that in the 
limit j4 ^ 0, one has 



lim E-a 



0. 
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It should be emphasized that although the energy per 
unit volume diverges as one approaches the initial cos- 
mic singularity, the maximum possible available energy 
at any given point is finite, and it goes to zero as t ^ 0. 

In order to study the behavior of the entropy per unit 
volume in the vicinity of the initial singularity, we can 
make an asymptotic expansion of the entropy for very 
small values of t. Defining the entropy density by s{t) = 
S'(t)/Vb, one obtains 
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+ ^E^Utf ] +0{l3{t)EraUtf)- 



While in the standard model f5{t) is given by ■Ait)- 
from our setting one can deduce that actually, must 
satisfy 
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The consequences of our result regarding the possible 
initial data for physical fields at the very early Universe, 
can be better understood if we recall from equation (|l|) 
that the entropy is given by 5 = — /cs Inp^; since this 

i 

expression implies that the only way in which the entropy 
of a system can go to zero is if all the probabilities pi go to 
zero except one, let us say poj which must have the unit 
value. This means that if we where to have two subsys- 
tems with identical statistical mechanic description and 
with zero entropy, then the physical state of both systems 
would be indistinguishable from the statistical mechanic 
point of view. This result explains way when we observe 
the cosmic background radiation in different directions 
on the sky that where causally disconnected at the time 
of decoupling, have such a similar description 

The attempts to understand the low entropy behavior 
at the early Universe has been the subject of study from 
different viewpoints. For example the standard presen- 
tations of the postulated inflationary Universes aims to 
explain the so called horizon problem, which is associated 
to the low entropy behavior ^ . 

According to our result we see that there is no prob- 
lem with the observed initial isotropic and homogeneous 
behavior of the Universe, regarding the existence of the 
initial cosmic singularity; since we have seen that such 
behavior is actually a consequence of the presence of the 
cosmic singularity. 
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which reduces to the usual behavior when /3(t) -Emax(0 
and /3(io) -E„iax(io) are big. Instead, when I3{t) Eraa^{t) 
is small, this expression shows a very different behav- 
ior for pit). From the fact that ii^max(i), for very small 
t, behaves as ^(i)^, one can deduce from equation ([T^ ) 
that Pit) behaves as Ait)^^ , in the vicinity of the initial 
singularity. 

Using both asymptotic behaviors; namely: ii^max(i) cc 
Ait)"^ and Pit) cc Ait)^'^, as t goes to zero, it is easily 
seen from equation ( |l2| ) that the entropy density goes to 
zero when one approaches the initial cosmic singularity; 
that is 
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